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ABSTRACT 

Neutrinos emitted during the collapse, bounce and subsequent explosion provide information about 
supernova dynamics. The neutrino spectra are determined by weak interactions with nuclei and nuclc- 
ons in the inner regions of the star, and thus the neutrino spectra are determined by the composition 
of matter. The composition of stellar matter at temperature ranging from T = 1 — 3 MeV and densi- 
ties ranging from 10 -5 to 0.1 times the saturation density is explored. We examine the single-nucleus 
approximation commonly used in describing dense matter in supernova simulations and show that, 
while the approximation is accurate for predicting the energy and pressure at most densities, it fails to 
predict the composition accurately. Wc find that as the temperature and density increase, the single 
nucleus approximation systematically overpredicts the mass number of nuclei that are actually present 
and underestimates the contribution from lighter nuclei which are present in significant amounts. 
Subject headings: dense matter - equation of state - supernovae: general 



1. INTRODUCTION 

Stars with masses larger than about 10 M Q end their 
lives in core-collapse supernovae. The initial collapse is 
initiated by the disappearance of the pressure support 
driven by (i) dissociation of nuclei and (ii) electron cap- 
ture on nuclei which make matter increasingly neutron 
rich. The electron captures emit neutrinos, which ini- 
tially escape, but are later trapped by the increasing den- 
sity in the core. The collapse is halted by the strong re- 
pulsion from high density matter and the core "bounces," 
driving an outward shock wave. The standard view is 
that this shock loses energy through continuing dissocia- 
tion of nuclei and neutrinos from the core finally restore 
energy to the shock driving an explosion. 

The evolution of the supernova is determined by the 
equation of state (EOS) and composition of matter at 
densities up to the nuclear saturation density, po ~ 
2.7 x 10 14 g/cm 3 and at temperatures between 1 and 
3 MeV. The collapse is driven by the tendency of m atter 
to have a low entropy per baryon (jBethe et al.lll979t ). the 
extent of electron ca ptures during col lapse is determined 
by the co mposition (|Hix et al.l 1200 3). and the neutrino 
spectrum (|Sumivoshi et al 120051 ) is determined by the 
nature of matter at the "neutrinosphere" , the surface of 
the last neutrino scattering, which is typically at 10 11 
g/cm 3 . 

At densities below p ~ 8 x 10 13 g/cm 3 , most matter 
resides in nuclei with A > 10 arranged in a Coulomb 
lattice, surrounded by a nucleon gas and embedded in 
a degenerate electron gas. Masses of isolated nuclei in 
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coexistence with the gas may exceed those of stable nu- 
clei but are limited to A < 1000 by the interplay of the 
Coulomb and symmetry contributions to their binding 
energies. Both nuclcons and nuclei together compose a 
nuclear (Fermi) liquid-gas phase mixture, whose Equa- 
tion of State (EOS) contributes significantly to the EOS 
for supernova simulations. 

A principal approach to computing the EOS for su- 
pernovae used in the past three decades has been 
the "Single-Nucleus Approximation" (SNA) where low- 
temperature matter is assumed to be composed of 
neutrons, protons, alpha part icles, and a single , 
representative, heavy nucleus (ILattimer et al.l 119851 ; 
ILattimer fc Swestvl 119911 ; IShen et all I1998D . The nu- 
clei vapor ize at temperatures t h at dep end sensitively on 
density. iBurrows fc Lattimerl (|1984l ) have found that 
the SNA can accurately predict thermodynamic func- 
tions for the matter, and microscopic calculations for 
the nuclear matter EOS have been performed in this 
approximatio n for densities rangin g from to satura- 
tion density (ILattimer et al.l 119851 ; ILattimer fc Swestvl 
119911 : IShen et alJ Il998| )" However, recent studies sug- 
gest that the SNA may be insufficient to completely 
describe th e dynamics of c ore-collapse supernovae. In 
particular, iHix et all (|2003h find that calculating elec- 
tron captures on a distribution of nuclei in nuclear sta- 
tistical equilibrium, in place of a single representative 
nucleus, increased th e potential for an explosion. Also, 
I Ar cones et al.l ()2008l ) find that the presence of light nu- 
clei is important for describing the neutrino spectra in 
the early post-explosion phase. 

A more realistic description of supernovae EOS would 
model the composition of matter by an ensemble of nuclei 
in near nuclear statistical equilibrium. In this work, we 
refer to this approach as the "Grand Canonical Approx- 
imation" (GCA). The GCA permits investigation of the 
shapes and widths of the nuclear mass and charge distri- 
butions and the population of nuclear excited states. In 
situations where the Wigner-Seitz approximation for a 
unit cell is valid and the widths of the heavy nuclear 
distributions can be neglected, GCA and SNA mod- 
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els should make simi l ar pre dictions. In a recent paper 
iBotvina fe Mishustinl (|2005l ) report on GCA calculations 
for supernovae and suggest that differences between their 
GCA and previous SNA calculations should be impor- 
tant. However, they did not directly compare SNA and 
GCA approximations to quantify the differences. 

To address this question, we ado pt a GCA model that 
is form ally equivalent to that of IBotvina fc Mishustinl 
(|2005f ). and construct an analagous SNA model for com- 
parisons to the GCA model. This corresponding SNA 
model employs the same mass formula and level densi- 
ties that are used in a corresponding GCA model. We 
compare predictions of the GCA model and its compan- 
ion SNA model for two different mass formulae, which 
differ by their treatment of the nuclear surface symmetry 
energy. Then, wc construct another SNA model (ISNA) 
which is m ore similar to the modern t abulated EOS con- 
structed in lLattimer fc Swestvl (|1991f ) and is frequently 
employed in supernovae simulations. The ISNA model 
includes interactions and the quantum statistics of the 
nucleons in the gas, a more complete description of the 
Coulomb energy, surface energy term which vanishes in 
the small proton number limit, a critical temperature 
which depends on the electron fraction, and effects asso- 
ciated with the presence of a neutron skin for nuclei with 
N > Z. 

From comparisons between the GCA and SNA mod- 
els, we find that the traditional SNA approach succeeds 
at most densities in describing the basic properties of 
the equation of state, but that it fails in describing the 
composition. We find that the differences between the 
predictions of GCA models and SNA models exceeds the 
differences between the predictions of the various SNA 
models. More specifically, we find that average quanti- 
ties, such as the mass fractions for neutrons, a particles 
and heavy nuclei, energy/baryon, and entropy, are sim- 
ilar at p < 0.01 and Y e = 0.4 in the GCA and SNA ap- 
proaches. The agreement is slightly worse for Y e = 0.2, 
and the agreement between GCA and SNA deteriorates 
even more for calculations that neglect the surface sym- 
metry energy. 

We find generally that the SNA and ISNA models agree 
quite well over most of the densities and temperatures 
considered; the largest differences are a few tenths of 
an MeV in the energy per baryon at higher densities 
and comparable differences in the entropy. The SNA 
and ISNA calculations, however, neglect completely the 
considerable widths of the mass and charge distributions 
predicted by the GCA calculations. At higher tempera- 
tures and densities, the SNA calculations systematically 
over-predict the size of the representative nucleus, and 
fail to treat the light nuclei that become more abundant 
in some regimes of density and temperature. These dif- 
ferences can have an impact on the electron capture and 
weak interaction rates that prevail in a supernova. 

In the following, we begin our comparisons by adopt- 
ing identical formula for the masses, level densities, and 
electron screening approximations for pairs of GCA and 
SNA calculations. We perform this comparison for two 
liquid drop models, LDM1, which has a surface symme- 
try energy term, and LDM2, which does not. Section 
2 describes these model assumptions. As these two cal- 
culations neglect the interactions and quantum statistics 
of particles in the gas, the ISNA model that contains 



these effects is described in Section 2 as well. Section 3 
describes and compares results obtained for the various 
GCA and SNA models. Section 4 summarizes this work 
and provides an outlook for future studies. 

2. THE MODELS 

Sections 2.1-2.3 describe the SNA and GCA models 
without the inclusion of the interactions and quantum 
statistics of particles in the gas. These two calculations 
can be compared to isolate the differences between the 
SNA and CGA models. Section 2.4 describes an im- 
proved SNA model (the ISNA model) that includes the 
interactions and quantum statistics of particles in the 
gas. This model allows tests of the relative importance 
of some simplifying assumptions used in the models de- 
scribed in Sections 2.1-2.3. 

2.1. Ground State Properties of Nuclei 

The binding energies of nuclei strongly influence the 
charge and mass distribu tions of nuclei wit h in hot ex- 
tende d nuclear systems (jSouza et al.l 120031 : iTan et all 
|2003|) . In this wor k we use a form for the liquid drop 
mass formula from (|Preston fc Bhadurilll975l ): 



B A . Z = C V A - C S A 2 / 3 - C c -^ + Sa^A- 1 ' 2 + C d ^ , 

where A and Z denote, respectively, the mass and atomic 
numbers, 



a 



1 - h 



A-2Z 
A 



(2) 



and i = v,s c orresponds t o volu me and surface, 
respectively. In ISouza et al.l (|2003f ). its parameters 
have been fitted to t he available experimental data 
([Audi fc Wapstr all 19951 ). The corresponding values are 
listed in Table [TJ and are labeled LDM1 (Liquid Drop 
Model). For completeness, the pairing term is also in- 
cluded in the expression above, although it is neglected 
in the calculations presented in this paper. It reads: 



Sa,z = 



-C, 



N and Z even 

A odd 

N and Z odd 



(3) 



where N = A — Z stands for the number of neutrons. Be- 
sides the standard terms, this parametrization includes a 
correction to the Coulomb energy due to the diffuseness 
of the surface, CdZ 2 /A, as well as a surface contribution 
to the symmetry energy. Both corrections are usually 
neglected in simple paramctrizations. While the latter 
contribution is most important for light nuclei, it also 
influences the masses of very heavy nuclei, as discussed 
below. 

The accuracy of this formula can be inferred from Fig- 
ure [TJ which shows the difference AB between the pre- 
dictions of equation (TTJ) and the empirical (measured) 
values. It is important to stress that AB corresponds 
to the difference between the total binding energies, i.e., 
it is not divided by the mass number. Since the above 
formula has no corrections associated with shell effects, 
pronounced discrepancies are observed near closed shells. 
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TABLE 1 

Parameters of the liquid drop mass formulas used in this work. All 

THE VALUES ARE GIVEN IN MeV. 
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Fig. 1. — Difference between the total binding energy predicted 
by the mass formula and the empirical values. For details, see the 
text. 

Shell effects for nuclei with A > 5 arc neglected since 
they do not strongly modify the qualitative shape of the 
nuclear distributions which we calculate. 

In order to investigate the influence of the Coulomb 
diffuseness correction and the surface symmetry energy, 
we have refitted the data while keeping the parameters 
Cd and ks equal to zero. In this limit, this mass formula 
has no surface symmetry energy or diffuseness correc- 
tion to the Coulomb interaction, which makes it have the 
same form that is used in many statistical calcu l ations 
dBotvina fc Mishustinl 120051. l2004t lAguiar et all 120061 ; 
iBhattacharvva fc Mekiianl 119991 ; iBondorf et all I1995D . 
The best fit parameters arc listed in Table [1] and labeled 
LDM2. The corresponding deviations AB from the em- 
pirical values are shown in the middle panel of Figure [TJ 
One can clearly see tendencies for LDM2 to over predict 
the empirical values for heavy masses and even stronger 
tendencies to under predict the empirical values for light 
masses. 

To aid in the discussion we define the nuclear symmetry 
energy, B sym : 



B K 



(A 



a s k s \ (A-2Z) 2 
-2Z) 2 



A 



(for LDM2). 



(4) 



The second line shows the simplification that occurs by 
setting k s = in the LMD2 paramctrization. 

The difference between the LDMl and LDM2 reflects 
the reduction of the symmetry energy in the nuclear sur- 



face and in the effective symmetry energy coefficient a sym 
in the LMD1 parametrization, which causes its masses to 
increase as a function of A, better following the trends of 
the measured masses. In the simpler LDM2 parametriza- 
tion, this coefficient is constant, which does not accu- 
rately follow the experimental trends and leads to the 
systematic deviations in both the light and heavy mass 
regions displayed in the figure. We have checked that the 
influence of the Cd coefficient does not account for the 
increase in AB observed in the heavy mass region. 

Experimental values are used for the very light nuclei 
(A < 5), whenever known, since both parameter sets give 
a very poor description of the binding energies of light 
nuclei. This procedure is not adopted for heavier nuclei 
as one has to consider nuclei far from the known mass re- 
gion. In this case, a careful extrapolation scheme should 
be devised so as not to introduce spurious effects in the 
calculated yields for mass re gions in which no experimen- 
tal in formation is available (jSouza et al.ll2003l : iTan et aTl 
l200l . Since the major thrust of this paper concerns the 
comparison of two approximations for calculating equi- 
librium distributions of nuclei, the use of a mass formula 
for all nuclei with A > 4 is sufficient and will ensure a 
smooth behavior for the predicted yields throughout the 
mass range considered below. 

A similar procedure is also adopted for the spin degen- 
eracy factors of nuclei. Empirical values are used only 
for A < 5. For heavier nuclei, we set the spin degeneracy 
factors to unity (vanishing nuclear spin). This approxi- 
mation is not that important because errors introduced 
by neglecting ground state spin degeneracies are much 
smaller than the uncertainties in the nuclear level den- 
sities at finite temperatures. These level densities are 
described below. 

2.2. The Grand Canonical Approximation (CCA) 

For a system composed of {k\, k%, ■ ■ ■ , A;m} different 
species in thermal equilibrium at temperature T, the 
grand-partition function is given by 



Zg = X) " ' Z c( T > V, h,k 2 , • • • , k M ) 

fei=0 k M =0 

e /3(/jifci+/i 2 *-'2H \-HKik M ) fg\ 

where f3 = 1/T, and the canonical partition function Zc 
is given by 



Zr 



9l V } A\' 2 /\% 



g M V f A% 2 /\% 



xe 



-0F(T,V,k u -,k M ) 



k M \ 



(6) 
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In the above expression, we ap- 

proximate the nucleu s translational e ffective mass by 
m A ,z = Am following iTan et all (|2003f ). We further ap- 
proximate m by the free nucleon mass m = 939 MeV/c 2 . 
The function F(T, V, k%, ■ ■ ■ , &;m) denotes the Hclmholtz 
free energy of the system, excluding contributions as- 
sociated with the translational motion, and gi denotes 
the spin degeneracy factor. The free volume Vf takes 
finite size effects into accou nt and is calculated as in 
iBotvina fc Mishustinl (|2004h . i.e., fl = V f /V = 1 - 
Pb/po, where V is the total volume, ps denotes the 
baryon density, and po represents normal (saturation) 
nuclear density. If F(T, V, ki, • ■ ■ , /cm) can be written as 

F = hh(T, V) + fc 2 / 2 (T, V) + --- + k M f M (T, V) , (7) 
the canonical partition function becomes 



Z C (T,V,k 



I,-'-, KM) 



M 

n 



[ 9l V f Af 2 /X 



which leads to 



M 



where 



Z G (T,V,{p*}) = Y[< 



-tr A 3/2 

q t - , 3 e 

A.rp 



(8) 
(9) 

(10) 



In this case, the average number of nuclei of species i in 
the volume V, Yi, can be easily evaluated 



oc kt oo ki oo k M 

Y - — V ^ V k — V 2m_ _ fl . 

ki=0 fci=0 km=0 

By assuming chemical equilibrium, the chemical po- 
tentials may be written as 

PA,z = PbA + p Q Z, (12) 
so that the number density of a given species (A, Z) be- 



comes 



V 



: ,( l i B A+fi Q Z-f A , z )/T 



(13) 



Thus, pb and pq are determined upon fixing the average 
number of baryons and by imposing charge neutrality, 
respectively: 



riB = y^^AnA.z 

A,Z 

n e = ^jZn A .z , 

A,Z 



(14) 
(15) 



where n P denotes the e lectron density. As in 
IBotvina fc Mishustinl (|2004l ), the sums include all nuclei 
with Z < A. We have checked that, for the temperatures 



and densities considered below, A = 1000 is a safe up- 
per cut in the sums above. We impose that the relative 
errors in equations (|14j|15p are smaller than 10~ 7 . 

In order to determine the chemical potentials, and sub- 
sequently, other relevant thermodynamical quantities, 
the Helmholtz free energy Ja.z must be specified. We 
include contributions to the internal excitation of nuclei 
associated with surface and bulk, besides the binding en- 
ergy of the nuclei and the Coulomb interaction among the 
particles: 



T- 



f A ,z=l3(T)A 2 / 3 A-B AiZ 

1/3 



Z z 
IV3 



n 



n 



where 



P(T) = A) 



T?-T 2 
T 2 + T 2 



5/4 



(16) 



(17) 



The values of the parameters (3q = 18.0 MeV, 
T c = 18.0 MeV, and e = 16.0 MeV corre- 
spond t o those usually adopted in statistical calcu- 
lations (lAguiar et all 120061: IBotvina fc Mishustinl 120041 : 
ISouza et al.ll2003l : ITan et al.ll2003l : iBondorf et al.lll995f) . 
The Coulomb energy, in the last term in equation (fl6|) . 
is calculated through the Wigner-Seitz approximation 
(jWigner fc Seitz| [l934V It should be noted that equa- 
tion (|16p differs from the usual expressions employed 
in statistical models because it takes into account the 
presence of the e lectron gas surrounding the nuclei 
(|Bavm et al.lll97ll ). In this respect, an important com- 
ment is in order. In principle, the baryon density in 
equation (|16|) should b e the average par ticle density in 
the Wigner-Seitz cell (|Bavm et al.l Il97lt ). However, if 
this is done, the Hclmholtz free energy cannot be cast 
in the form of equation ([7]) and, then, the simple for- 
mulae above would n o lon ger be valid. Thus, as in 
IBotvina fc Mishustinl (|2004h . we approximate the aver- 
age density by ns in order to keep the simple relations 
above. Since we want to maintain the GCA and the 
SNA close as possible, we adopt this approximation in 
SNA model as well. 

Finally, nuclei with A < 4 are treated as point parti- 
cles, with no internal degrees of freedom. Therefore, the 
temperature dependent terms in equation (|16[) are set to 
zero for those nuclei. 

2.3. The Single Nucleus Approximation (SNA) 

Since the determination of the chemical potentials 
above and the subsequent computation of the rele- 
vant thermodynamical quantities is too time consum- 
ing to be used in many practical ast rophysical calcula- 
tions , the GCA has been s impli fi ed (lEpstein fc Arnettl 
19751: iBurrows fc Lattimen 119841 : lLattimer et all 119851 : 
Lattimer fc Swestvlll99lh . More specifically, instead of 
considering all possible species, only a few types are al- 
lowed: neutrons, protons, alpha particles and a single 
species of heavy nucleus. The latter should, to some ex- 
tent, represent the nucleus distribution at the heavy mass 
region. The other three nuclei are intended to consider 
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the main contribution to the nucleus distribution in the 
light mass region. 

The constraints associated with baryon number and 
charge neutrality simplify to 



n B = n n + n p + 4n a + A h n h 



(18) 



and 



n p + 2n a + Z h n h 



(19) 

where n n , n p , n a , and rih respectively denote the num- 
ber density of neutrons, protons, alpha particles, and the 
heavy nucleus of mass and atomic numbers Ah and Zh- 
In this way, the most probable configuration is found by 
minimizing the total free energy of the system T ', subject 
to the above constraints 



d 

— VXi [n B - n n -n p — 4n Q - A h n h ] 

- VX 2 [tie -n v - 2n a - Z h n h }} = . (20) 

In this expression, A; stands for the Lagrange-multiplicrs 
and Xi denotes Ah, Zh, n n , n p , n a , and nh- The to- 
tal Helmholtz free energy possesses the same ingredients 
used in the GCA: 



+ (f a + f a rans )n a + (fh + ffrn^h] . (21) 

The contribution due to the translational motion is given 

by 



-T 



log 



and, for large values of k%, one may write 



ft 



= -T 



log [ I^La^ 2 



1 



(22) 



(23) 



since log(fcj!) m k t \og(k t ) - k { . 

The derivatives associated with n n and n p allow one 
to easily eliminate the Lagrange multipliers: 



Ai =Tlog 



A 3 77 



(24) 



A 2 = Tloe 



3 (n B 
2 



1/3 



1 fn B 

2 I no 



(25) 

The remaining derivatives, with respect to Ah, Zh, n a , 
and n /j, lead to non-linear equations, which must be 
solved numerically Due to the logarithmic factors en- 
tering in the formulae, it is convenient to use 



^ J 4 3/2 e - (/, " fI,,/T 



(26) 



so that the equations are solved for the chemical poten- 
tials of each species. We require the same precision on 
the constraints (|18H19[) as in the grand canonical calcu- 
lations. 



2.4. A Single Nucleus Approximation with 
Interactions (ISNA) 

This improved s ingle nucleus model (ISNA) uses the 
mass formula from Stcincr ( 2008) , adapted slightly to in- 
clude finite temperature, alpha particles, and the effects 
of the neutron skin. In this model, the binding energy of 
a nucleus at temperature T , with proton number Z and 
mass number A, in the presence of an external proton 
gas of density n PjOUt , is given by 

B(Z,A,n n ,n p ,n PtOUt ,T) = B hulk (n n ,n p ,T) 

+aB{n n ,n p ,T) 



+C e c 



oulomb 



\ 77 2 

(Z,A,n p , out )/(An). 



Here, n n , n p and 77 = 



(27) 

denote the average in- 
ternal neutron, proton and baryon number densities, re- 
spectively. 

The binding free energy of bulk matter, -Bbuik is given 

by 

-tfbuik = [e(n n , n p ,l) - n n m n - n p m p 



-Ts( n 
-Ts(n n ,0,T)] 



Ar 



skin 



[e(n n ,0,T) - n n m T 



(28) 



where m n and 777 p are the neutron and proton masses, 
e(n n ,n p ,T) is the energy density of homogeneous mat- 
ter evaluated at the given neutron and proton density 
and temperature, and s(n n ,n p ,T) is the entropy den- 
sity. The energy density of homogeneous matter is de- 
scribe d using the Akmal , Pandharipande, and Ravenhall 
EOS (jAkmal et al.lll998l . APR). Finite temperature cor- 
rections to the APR EOS are computed assuming that 
there are no finite temperature corrections to the poten- 
tial energy contribut ions similar to the approach used in 
iPrakash et alJ (|1997l ). The part proportional to A s km is 
a correction to the bulk energy density for the neutron 
skin. 

The average baryon density in equation (|28[) is deter- 
mined from 

n = n n + n p = n s + ndl 2 (29) 
where 1 = 1 — 2Z/A. The parameter n s is analogous to 
the saturation density of nuclear matter and is expected 
to be near 0.16 fm -3 . The parameter rid (which is typ- 
ically negative) subsumes the decrease in the saturation 
density with the isospin asymmetry. 

The individual average neutron and proton number 
densities are given by 

n„=n(l + 5)/2 

n p = n(l-6)/2 (30) 
and the density asymmetry S = l — 2n p /(n n +n p ) is given 
by 5 = £1 where C is a constant parameter determined 
by the fit to the experimental nuclear masses. 
The parameter A s ki n is chosen so that 



Air 

A - A^Hn = -^-R p (n n + n p 



(31) 



The neutron and proton radii, R n and R p are fixed by 

47T 



A-Z-- 



-R n n n 



47T 3 

Z = —R p n p . 



(32) 
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In summary, for the bulk part of the nuclear mass for- 
mula, using Z and A one can compute the average neu- 
tron and proton densities using the relations above, and 
the radii from equations (|32[) , then compute iV s kin from 
equation (|31|) . and use the equation of state for homoge- 
neous matter to compute the bulk part in equation (|28| 
above. Note that the symmetry energy contribution to 
the nuclear mass is automatically counted above as part 
of the bulk contribution. 

The surface energy density as a function of the "surface 
tension" , a is 

^surface = ~~ (33) 
K 

For the mass formula, we need the surface energy per 
baryon, which for T = and / = is given by 



B(T = 0,1 = 0) surface M 



da 

TlR' 



Where the radius, R is determined from 
4ttR 3 



-n = A 



so that 



B(T = 0,1 = 0) surfacc /A 



3cr / 47m \ 
~ \~3a) 



1/3 



(34) 



(35) 



(36) 



In general, the surface energy should be modified to en- 
sure that the surface energy vanishes in th e limit S — > 1 as 
it mus t. To address this issue, we follow lLattimer et al.l 
(|1985f ) and approximate B{n n ,n p ,T) by: 



B(n n , n p , T) = r\(n n ,n p , T) 



16 + b 



[l/x^ + b+l/(l-x) 3 } 
xB(T = 0,/ = 0) surfacc , (37) 

where x = n p /n and 6 is a free parameter. 

For T > 0, one must consider the reduction in the sur- 
face tension of the nucleus due to the interactions with 
the surrounding gas. For the finite te mperature correc- 
tion, we follow (|Lattimer et al.lll985l ) and approximate 
r]{n n ,n p ,T) by: 



i](n n ,n p ,T) 



(i-r 2 /rg) 

(1 + aT 2 /^) 



5/4 



(38) 



This expression is essentially identical to the SNA 
model except for the presence of the factor a = 
0.935 - 5.1(0.5 - xf - 1.1(0.5 - x) 4 which is ob- 
tained from a fit to the results of iRavonhal l et al.l 
(|1972h that allows extrapolations to neutron rich sys- 
tems. The critic al temperature is also isospin depen- 
dent. Following lLattimer et alj (|1985l ) . we extrapo- 
late Tc to neutron rich matter using Tc = Tq{x = 
1/2)^/1 - 3.313(0.5 - x) 2 - 7.362(0.5 - x) A . We also 
take T c (x = 1/2) = 20.085 MeV as was done by 

lLattimer etTaTl (|1985l ). 

The Coulomb energy density (|Ravenhall et al.l Il983f ) 



is 



^Coulomb 



:7r (n p - n Pi0 ut) e R 



1/3 



x P 



(39) 



where e 2 is the usual Coulomb coupling ~ Hc/137 and 
Xp = Rp/Rfi is the volume fraction of matter occupied 
by the proton sphere. The Coulomb contribution is mul- 
tiplied in equation (|27p by a final parameter, C, which 
takes into account the fact that the proton density dis- 
tribution has its own surface. The parameter values are 
given in Table CD and labeled LDM3. 

In order to determine the composition and properties 
of matter in the supernova, we minimize the free energy 
at a fixed density as a function of the proton number 
and atomic number of nuclei, and the number density of 
dripped neutrons, n n ^ r ip- The free energy of this matter 
is given by 

f {Z , A, Tl n ^ ou i, Tl p ^ ou t, Tl a ou i, T^j 

nB(Z, A, n„ iOUt , n PiOUt , T)/A 
+(1 - X) [/ (»Vout, %,out, T) 

+ f(n a ,out, T)} + /classical ("Nuclei, T) + f c \(n e , T$40) 

where "Nuclei = X n /A, x is the volume fraction of matter 
occupied by nuclei, and /classical ("Nuclei, T) is the classi- 
cal free energy from the ideal gas of nuclei, and f \(n e ) 
is the free energy of the electrons. For the purposes of 
comparisons with the SNA and CCA models, we do not 
include the free energy of the electrons in the ISNA cal- 
culations presented in this paper. The free energy of the 
dripped nucleons, /(" ra>0 ut> "p,outj T), is computed with 
the APR EOS in the same way as the bulk contribution 
to the nuclear mass formula referred to above. 

The constraints of baryon and charge conservation are 
implemented with 

n B = Xn n + X P n p + (1 - x) "n,out + (1 - Xp) n PlOU t 
+4(1 -x) n a 

Y e n B = Xp n p + (1 - Xp) "p.out + 2 (1 - x) "a (41) 

One of the two constraints is used to fix x, and the other 
can be used to constrain one of the parameters to the 
free energy, e.g. n Pi out- 

3. RESULTS AND DISCUSSION 

We utilize the models described above to compute 
the composition and thermodynamic functions for mat- 
ter for the densities 10~ 5 < "b/"o < 10 _1 , temper- 
atures 1 MeV < T < 3 MeV, and electron fractions, 
0.2 < Y e < OA relevant for core-collapse supernovae, see 
e.g. Uanka et all (|2007f) . 

The isobar number density 



n A 



z 



TiA.Z 



(42) 



is calculated through equation (Tj"3")) using the LDM1, and 
the results are displayed in Figure [5] for different baryon 
densities at T = 1.0 MeV and for two representative val- 
ues of the electron fraction used throughout this work. 
The lines correspond to the distribution predicted by the 
GCA, whereas the arrows indicate the value of Ah in 
SNA. The results reveal that A^ is systematically larger 
than the mass number A max at which the maximum of 
the isobar number density occurs, in the heavy mass re- 
gion. This may be explained by noticing that heavy nu- 
clei, many of them heavier than A max , contribute sig- 
nificantly to the sums in equations (|14[) and (TT5)) due 
to the smooth behavior of ha z- These contributions 
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500 



Mass number 



Fig. 2. — Isobar number density for different values of ng at 
T = 1.0 MeV, using the LDM1 mass model. The lines represent 
the results obtained with the GCA, whereas the arrows point to 
the values of A h predicted by SNA at the same values of iij used in 
the GCA. The baryon densities associated with the arrows increase 
from left to right. 

also include lighter nuclei such as the hydrogen isotopes 
(d,t), helium isotopes ( 3 He, 6 He) and intermediate mass 
nuclei with 3 < Z < 20. These other contributions, 
combined with the overall mass and charge conservation 
constraints, lead to A max values that are lower than Ah- 
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Fig. 3. — Mass fractions associated with nuclei with Z = 0, 
Z = 2, 3 < Z < 20 (intermediate), and Z > 20 (heavy), obtained 
employing the LDM1. The solid, dashed and dotted curves denote 
the predictions of GCA, SNA and ISNA calculations. Upper, mid- 
dle and lower panels provide the mass fractions for neutrons, alpha 
particles and heavy nuclei. The dashed-dotted curves in the mid- 
dle left and middle right panels denote the mass fractions in the 
GCA calculations due to nuclei with 2 < Z < 20, excluding alpha 
particles. The left panels give results for Y e = 0.2, and the right 
panels give results for Y e = 0.4. All panels are for T = 1.0 MeV. 



It may also be noticed that Ha shows a broad distribu- 
tion for all baryon densities displayed in the figure. The 
width becomes narrower as decreases, but it remains 



finite, for the reasons discussed above. The position of 
the peak of the distribution also shifts to lower values 
since dilute configurations favor partitions with higher 
numbers of free nuclcons and light particles, which com- 
pete for the available charge and mass (equations (|14p 
and (fl"5j) ) against heavy nuclei. 

This competition is illustrated in Figure [31 which 
shows the mass fraction associated with different par- 
ticles. These results clearly reveal that more mass is 
contained in neutrons and alpha particles as one goes 
towards lower densities. The contribution from Z = 1 
is not presented since it is relatively small for the neu- 
tron rich systems we consider (even though it is included 
in our mo dels). Consistent with the conclusions of ear- 
lier work (lEpstein fc ArnettJ 119751 ; iBurrows fc Lattimerl 
Il984t lLattimer fc Swestvlll99lh . the mass fractions of al- 
pha particles, heavy nuclei, and neutrons predicted by 
the GCA and SNA calculations are very similar, despite 
the non-negligible widths of the mass and charge distri- 
butions for the CGA calculations. 

Similarly, the interactions between the gas particles 
and between the gas and the heavy nucleus included in 
the ISNA model make comparatively small differences 
to the mass fractions in Figure [3] The main differences 
appear to be that the ISNA model predicts somewhat 
smaller values for alpha particle and larger values for 
the neutron mass fractions compared to the SNA model. 
This partly reflects small differences between the mass 
formulae used in the SNA and ISNA models. In general, 
both models predict similar values for the the mass of 
the heavy nucleus at Y e = 0.4, but the heavy nucleus is 
somewhat lighter for the ISNA model at Y e = 0.2. 

The qualitative dependence of A max and Ah on the 
electron fraction may be understood in terms of the sym- 
metry energy. At small values of Y e , neutrons are more 
abundant than protons and consequently, the neutron 
chemical potential exceeds the proton chemical potential. 
As a consequence, the number of free neutrons and the 
neutron mass fraction increase, as can be noticed from 
Figure [3l and nuclei also tend to have more neutrons 
than usual. Owing to the symmetry energy, an increase 
in the neutron number increases the binding energy of 
protons within a nucleus, the tendency for nuclei to have 
larger neutron number N is accompanied by the distri- 
bution of heavy nuclei shifting to larger A and Z values. 
The increased neutron mass fraction at Y e = 0.2 can only 
be achieved if the heavy mass fraction correspondingly 
decreases. For this reason, the mass fraction associated 
with heavy nuclei is smaller for Y e = 0.2 compared with 
those obtained at Y e = 0.4. 

At higher densities, the decrease in the free volume 
and corresponding increase in heavy nucleus mass means 
that the additional neutrons at Y e = 0.2 will be largely 
contained in the heavy nucleus. Thus the isotope com- 
position of the heavy nucleus should be more sensitive 
to the baryon density when the matter is appreciably 
asymmetric. This aspect is illustrated in Figure IH which 
shows the average value of Z/A obtained with both mod- 
els. As expected, it is fairly independent of ub for nearly 
symmetric matter, whereas it is quite sensitive to it at 
small values of Y e . The predictions of the SNA model 
follow those made by the GCA fairly closely, at least 
at low densities. The SNA model tends to give more 
neutron rich heavy nuclei, partly because Ah is system- 
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Fig. 4. — Average Z/A of heavy nuclei predicted by the GCA, 
SNA and ISNA models, using the LDM1, are shown by the solid, 
dashed and dotted lines, respectively. The dashcd-dotted lines 
show the average Z/A values for lighter nuclei with 3 < Z < 20 
that lie below the heavy nucleus charge distribution. 

atically larger than A max and the isobar with the highest 
binding energy in nuclear mass formulae is increasingly 
shifted towards the neutron rich side with increasing nu- 
clear mass. Some of the extra neutrons in the GCA cal- 
culations are also contained in lighter intermediate mass 
nuclei with 3 < Z < 20 that lie below the heavy nucleus 
charge distribution. 

The collapse of a supernova is accompanied by large 
changes in many properties of the system, such as the en- 
tropy , pressure, temperature, density, and energy (|Bethd 
Il990h . some of which have an impact on the supernova 
dynamics. We therefore have evaluated the entropies and 
energies predicted by the GCA, SNA and ISNA mod- 
els. The entropy may be obtained from these approaches 
through the standard thermodynamic relation 



S = log Z G + T 



9 1 7 
— - log Zn 

dT 



(43) 



Applying the relationships in the preceding section, S 
becomes 



S = V^n AiZ 

A,Z 

where 



5 , 1 , df A ,z 
2 + j,UA,z - VA,z) q^t 



TT 2 



I C 



2\ V4 



9 f A > z = Ka A 2/3_ 
dT ' ° (T 2 + T 2 ) 2 \T 2 +T 2 

The total energy of the system is given by 
E = VJ2 nA,z»A,z + TS- PV. 

A,Z 

The pressure P, which may be obtained from 



(44) 



2A 

— T . 

(45) 



(46) 



P = T 



w log Zo 



E 

A,Z 



n A ,z 



in z 2 

2TA~^ 



CrU 



where 



(n s /n ) 1/3 + (ns/no) 



(47) 



(48) 



The second term in equation (|47|) is due to the Coulomb 
interaction between the electron gas and the nuclei and, 
therefore, is always negative for ns < no- 




Fig. 5. — Entropy and energy per baryon for different values of 
the baryon density, obtained using the LDM1. The solid, dashed 
and dotted lines denote the CGA, SNA and ISNA calculations 
respectively. The left panels are for Y e = 0.2 and the right panes 
are for Y e = 0.4. 



The entropy and energy predicted by the GCA and 
SNA models are shown by the solid and dashed lines in 
Figure \E\ for the same baryon densities and temperature 
used above. We do not show the nuclear pressure since, 
at these d ensities, it is much smaller than that of the elec- 
tron gas (|Bethejll99fi ICooperstein1ll985D and, therefore, 
does not play a significant role. Although a fairly rea- 
sonable agreement is found for low densities ub < 10~ 3 , 
differences between GCA and SNA models can be found 
at higher densities where the GCA mass distributions 
become rather broad. The deviations between the SNA 
and ISNA results for the energy and entropy are rela- 
tively small, only a few tenths of an MeV for Y e = 0.4. 
However, there are larger deviations at Y e = 0.2 where 
the nuclei are extremely neutron rich. 

The differences at higher densities and low Y e are 
not surprising: this is the region where one expects 
the results to depend sensitively on the details of the 
mass formul a and informa tion about the nuclear symme- 
try energy (|Steinerl l2008f ). However, the differences at 
low density have a different origin: near Y e = 0.2 and 
ns/no = 10 -5 , the free energy is very flat in the direc- 
tion of the parameter space which determines the num- 
ber of nuclei relative to the number of alpha particles. 
The reason for this is that the free energy can be mini- 
mized in one of two ways: (i) the system can choose to 
create more entropy (and thus decrease the free energy) 
by making a lot of alpha particles from nuclei, or (ii) 
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the system can choose to create nuclei from alpha parti- 
cles because the extra binding created by nuclei will also 
decrease the free energy. Within the ISNA model one 
can vary the number density of nuclei by 70 %, and this 
modification only changes the free energy by 0.1 %, while 
changing the entropy by 10 %. Because the free energy 
is so flat, the entropy in this region is quite sensitive to 
the nuclear mass formula. Fortunately this low-density, 
low-y e region is not actually probed frequently in actual 
supernova simulations. 




Mass number 



Fig. 6. — Same quantities as in Figure[2] In this figure, however, 
the results are obtained using the LDM2 mass model. 

In order to verify the extent to which our conclusions 
depend on treatment of the nuclear symmetry energy in 
the liquid drop formula used in these statistical treat- 
ments, we have also carried out GCA and SNA calcula- 
tions using the LDM2 formula presented in Sect. 12. ll The 
corresponding isobar number densities are shown in Fig- 
ure [S] Qualitatively, the calculations in Figure [S] resem- 
ble those in Figure [2] obtained using the LDM1 formula. 
Quantitatively, there arc differences, however. Most 
notably, the discrepancies between the GCA and SNA 
model predictions are now much larger for LDM2 than 
for LDM1, particularly for Y e = 0.2, at ns/no = 0.1. 
This is due to the lack of a surface symmetry energy 
term in the LDM2 mass model. The resulting under pre- 
diction of the binding energies for light masses in the 
LDM2 formulae shifts mass of the maximum A max for 
the GCA distribution towards larger values, where the 
difference between A max and Ah values is typically much 
larger. Nevertheless, the mass fractions predicted using 
the LDM2 mass model (not shown) arc similar to those 
shown in Figurc[31 except that intermediate mass nuclei 
(3 < Z < 20) give a somewhat larger contribution for 
the LDM2 to the mass fraction at the lowest density. 

The discrepancies between the predictions of the treat- 
ments are more pronounced at the highest baryon den- 
sity. To illustrate the dependence on temperature there, 
we compare the two calculations at ns/no = 0.1 for tem- 
peratures ranging from 1 to 3 MeV, values that are in the 
range relevant to supernova studies. The isobar number 
densities, obtained using the LDM1, are shown in Fig- 
ure [7] Reflecting the larger phase-space accessible to the 
system at higher temperatures, the GCA isobar distri- 
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Fig. 7. — Same as Figure [2] except that the baryon density is 
fixed and the temperature is varied from 1 to 3 MeV. 



butions become broader, the value of A max decreases, 
and a larger fraction of the mass is in the form of nuclc- 
ons or light nuclei at higher temperatures. The trend of 
increasing nuclcon yields with temperature is also exhib- 
ited by the SNA model, but the discrepancies between 
the GCA and SNA approaches become larger as T and 
the widths of the mass distributions increase. In spite of 
the larger width of the isobar distributions, however, the 
mass fractions obtained with GCA and the SNA model 
(not shown) agree to within 15% for Y e = 0.4 and to 
within about 25% for Y e = 0.2. The intermediate mass 
nuclei have a slightly larger mass fraction, but not enough 
to significantly change the mass fraction associated with 
heavy nuclei. 




Mass number 



Fig. 8.— Same as Figure [3 except that the T = 3 MeV. The 
inset shows values for calculated within the SNA model as a 
function of the density. 

Figure [8] shows the complementary comparison be- 
tween the GCA and SNA calculations using the LDM1 as 
a function of density at T — 3 MeV. We do not compare 
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the two calculations below n/no = 4x 10~ 3 because at 
this temperature the mass Ah of the most probable heavy 
nucleus drops to zero below this density. The inset in the 
figure shows the detailed dependence of Ah predicted by 
the SNA model. Even though, some heavy nuclei are 
predicted by the GCA model at n/n < 4 x 10 -3 , the 
mass fraction of Z > 20/ A > 50 at such densities is small 
and of the order of 3 — 6 x 10~ 2 . Both Figure and Fig- 
ure [8] show that the discrepancies between the GCA and 
SNA calculations increase with density over the range of 
densities investigated in this work. 



at densities of n/no ~ 0.01 — 0.1 lead to non-negligible 
differences between the entropies and energies predicted 
by the GCA and SNA approaches. The ISNA results are 
similar, again showing larger deviations from SNA calcu- 
lations in the neutron rich case at Y e = 0.2. The density 
range over which results arc available is slightly smaller 
because the disappearance of the heavy nucleus occurs 
as a higher density in the ISNA model. The transition 
from heterogeneous to homogeneous matter, especially 
in neutron-rich matter, is very model dependent. 
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Fig. 9. — Same quantities as plotted in Figure [5] except that the 
baryon density is fixed at p/po = 0.1 and the temperature is varied 
from 1 to 3 MeV. 
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Fig. 10. — Same as Figure [5] except that the temperature is fixed 
at T = 3 MeV and the baryon density is varied. 



In general, the agreement between the GCA and SNA 
approaches is somewhat better at Y e = 0.4 than at 
Y e = 0.2 where the mass distributions are wider. This 
conclusion also holds for the thermodynamic variables. 
The dependence of the entropy and energy per baryon on 
the temperature is shown at a fixed density of n = O.lno 
in Figure [5] Clear differences between the GCA and 
SNA calculations are predicted for Y e = 0.2. There the 
entropy difference (AS/A as 0.2) remains roughly con- 
stant as a function of temperature. The larger mass in 
the SNA leads to a lower energy baryon in the SNA than 
in the GCA calculation. The energy difference reaches 
its largest value (AE/A w 0.15 MeV) at T = 1 MeV. 
This difference decreases with temperature to a negli- 
gible value at T = 3 MeV. In contrast, the differences 
between GCA and SNA calculations are comparatively 
small for Y e = 0.4, reflecting the smaller widths of the 
heavy mass distributions at Y e = 0.4. 

The complementary density dependencies of the en- 
tropies and energies for GCA and SNA calculations at 
T = 3 MeV arc shown in Figure [TO] For both Y e = 0.2 
and Y e = 0.4, the largest differences between the ener- 
gies of SNA and are actually observed at n/n « 0.01. 
The entropy differences generally increase with density 
for both Y e = 0.2 and Y e = 0.4. Both entropy and en- 
ergy differences are larger for Y e = 0.2 than for Y e = 0.4, 
reflecting the larger widths of the heavy nucleus mass dis- 
tributions at Y e = 0.2. Both Figs. 151 and ITD1 demonstrate 
that the wide mass distributions of the GCA calculations 



Finally, to provide one more example of the importance 
of the mass formula, we show, in Figure 1111 the mass 
distributions predicted by the GCA for the LMD2 mass 
formula in comparison to the corresponding SNA predic- 
tions. We find that there are larger discrepancies between 
the GCA and SNA predicted for the LDM2 parametriza- 
tion than for the LMD1 paramctrization. As the mass 
fractions and the thermodynamic variables obtained in 
this case are very similar to those obtained for the LDM1 
paramctrization, we do not show them in the interest of 
brevity. Comparing the predictions in two different pa- 
rameterizations of the nuclear symmetry in Figure 1111 
we find that the changes in the mass distributions due to 
inclusion of the surface symmetry energy correction are 
less significant than the predicted differences between the 
GCA and SNA calculations that use the same liquid drop 
mass model (and same symmetry energy). 

In principle, nuclei must have a surface symmetry en- 
ergy term. The magnitude of the surface symmetry 
energy coefficient, however, is not well constra ined by 
the measured nuclear masses (jDanielewiczl I2003T ) . Even 
though the differences between the GCA and SNA cal- 
culations exceed the differences between the LDM1 and 
LDM2 calculations, parameterizations for the density de- 
pendence of the symmetry energy may be chosen tha t 
predict much larger effects (jBotvina fc Mishustinll2005T ). 
Ideally, a general investigation of the role of the symme- 
try energy in the supernovae EOS should include both 
an assessment of the sensitivity of the mass distributions 
to the choice of the symmetry energy as well as an as- 
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sessment of whether the chosen f orm is consistent w ith 
known experimental information (|Souza et al.ll2008h . A 
systematic investigation of the effect of the symmetry 
energy in relationship to the current experimental infor- 
mation is in preparation. 
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Fig. 11. — Same as Figure [7] except that the LDM2 is used 
the calculations. 



4. CONCLUDING REMARKS 

This work has focused on the differences between 
the Grand Canonical Approximation and simplified 
Single Nucleus Approximation utilized in most super- 
nova calculations. Consist ent with the earlier work in 
iBurrows fc Lattimerl (|1984f ), we find that the thermo- 
dynamic functions are relatively well described by the 
SNA at densities below O.Olpo- At higher densities, dif- 
ferences emerge that can be traced to the large widths 
of the GCA mass distributions. We have compared the 
simple SNA model to a more elaborate ISNA model that 
considers effects that become relevant at higher densities, 
such as the interactions between gas particles. We find 
differences between the SNA and ISNA models of these 
average quantities that are comparable to the differences 
between the GCA and SNA calculations. This indicates 
the importance of taking such effects into account accu- 
rately. Comparisons of the GCA, SNA and ISNA models 
suggest that the error made by assuming a single nucleus 
in the energies and entropies per baryon may be compa- 
rable to or smaller than the other uncertainties in the 



problem, such as the nuclear mass formulae and the in- 
clusion of quantum statistics and interactions between 
gas particles. 

Large differences can be found between the GCA pre- 
dictions for the the masses of nuclei and those of the 
SNA and ISNA calculations. Single nucleus approxima- 
tions, in general, tend to predict much larger heavy nu- 
cleus masses than does the GCA. Considering first the 
differences between GCA and SNA calculations, we find 
that the SNA model systematically ovcrprcdicts the mass 
number of nuclei that are actually present in the more 
accurate GCA model because it neglects the contribution 
from lighter nuclei and the large width of the mass dis- 
tribution. This discrepancy grows with increasing tem- 
perature. Such a large difference in nuclear distributions 
is significant. Even though some regimes (particularly 
higher densities and lower temperatures) are affected by 
quantum statistics, our calculations indicate that the ef- 
fects of assuming a distribution of nuclei rather than 
a sole representative nucleus are much larger and suf- 
ficiently robust. This suggests that first explorations of 
the consequences on supernova simulations of these dif- 
ferences using classical models of non-interacting nuclei 
may be helpful and should be performed. 

These results has specific implications for the equation 
of state tables and the weak interaction physics inputs 
used as inputs for supernova simulations. Our work sug- 
gests that tables utilizing the SNA can be updated to 
include nuclear distributions by applying an appropri- 
ate GCA model. We are exploring such issues with the 
aim o f determining their impact on the weak interaction 
rates (|Janka et al.ll2007t ) in supernova simulations. The 
use of a GCA model will also allow future investigations 
to more accurately determine the neutrino spectra be- 
cause it directly handles light nuclei, the imp ortance of 
which was emphasized in lArcones et al.1 (|2008| ) . 
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